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We develop an alternative formulation of the theory of solitons, polarons and multipolarons 
in quasi-one-dimensional degenerate and non-degenerate conducting polymers, starting from the 
continuum Hamiltonian introduced by Brazovskii and Kirova. Based on a convenient real-space 
representation of the electron Green function in one dimension, we present a simple method of cal- 
culating the Green function and the density of states in the presence of a single soliton or polaron 
defect, using which we derive exact expressions for the soliton, polaron and multipolaron excitation 
energies and the self-consistent gap functions for an arbitrary value of the electron-phonon coupling 
constant. We apply our results to cis-polyacetylene. 



On 
OS 



Q 



Os 



o 
o 



X 



I. INTRODUCTION 



> 

A wide variety of fascinating physical phenomena occur in quasi-one-dimensional materials with a Peierls-Frolich 
ground state [0J^]- In particular, the proposal ||-|7j that non-linear excitations such as solitons and polarons play a 
crucial role in the electronic properties of conjugated conducting polymers and other related materials have attracted 
great interest over the past two decades |2|,|8| O. Experimental and theoretical efforts to confirm the existence of 
these non-linear excitations and to clarify their properties are being continued to the present day |]l^-]l6| . 

In a large number of theoretical studies of solitons and polarons, the electron-electron interaction is ignored |T^ ] 
and the lattice motion is treated classically. A simple lattice model of non-interacting electrons in one dimension 
coupled to phonons was introduced by Su, Schrieffer and Heeger Q and has been applied mainly to numerical studies 
of irans-polyacetylene [(CH) X ], which is a representative conjugated polymer with two degenerate ground states. 
A continuum version of this model derived by Takayama, Lin-Liu and Maki [18] admits exact soliton and polaron 
solutions and has been the starting point in a number of analytical studies [Q, 19 j. In order to treat non-degenerate 
materials with a unique ground state and a higher-energy metastable state such as czs-polyacetylene, Brazovskii and 
Kirova derived a generalized version of the continuum model, which contains a new parameter, Ag, representing the 
extrinsic gap in the electronic spectrum that exists even in the absence of the Peierls distortion |6|. When A e = 0, 
Brazovskii and Kirova's model reduces to Takayama, Lin-Liu and Maki's model. If A e ^ 0, it does not admit the 
soliton solution, but allows the so-called multipolaron solutions in addition to the polaron solution. 

In this paper, we reformulate the theory of solitons, polarons and multipolarons, starting from the continuum 
Hamiltonian introduced by Brazovskii and Kirova, which can describe both degenerate and non-degenerate conduct- 
ing polymers depending on the value of A e . We develop an efficient method for calculating the Green functions 
associated with the non-linear excitations, based on a convenient real-space representation of the Green function in 
one dimension. Using the density of states expression obtained from the Green function, we compute the soliton, 
polaron and multipolaron excitation energies and the related self-consistent gap functions exactly for an arbitrary value 
of the electron-phonon coupling constant. These results are usually presented in the weak electron-phonon coupling 
limit in the literature and, to the best of our knowledge, have never been written down explicitly before. Our method, 
which involves only concepts familiar in condensed matter physics, is simple and transparent and has the advantage 
of being readily generalizable to more difficult problems such as the influence of disorder on solitons and polarons. 
In the next section, we introduce the Hamiltonian and the soliton wave function. The polaron case will be discussed 



in Appendix A. In Sec. [II and Appendices A and B, we describe our method for calculating the Green function and 
the density of states. In Sec. [TV], we compute the soliton, polaron and multipolaron excitation energies and the 
self-consistent gap functions and apply the results to cis-polyacetylene. In Sec. we conclude the paper with some 
remarks. 
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II. HAMILTONIAN AND WAVE FUNCTIONS 



We consider the continuum Hamiltonian of quasi-one-dimcnsional conducting polymers and related materials first 
introduced by Brazovskii and Kirova f(||2(| : 

dx A*(x), 



1 



(1) 



where vf is the Fermi velocity and A is the dimensionless electron-phonon coupling constant (For a precise definition 
of A in terms of the parameters of the lattice model, refer to a^s (i — 1,2,3) are Pauli matrices. The (real- 

valued) gap function is written as A(x) = Ai(x) + A e , where A^a;) is the intrinsic gap function, which is sensitive to 
electron-phonon coupling and varies slowly on the scale of a lattice spacing, and A e is a constant extrinsic component. 
The time derivative of Ai(x) is assumed to be sufficiently small that the lattice kinetic energy can be ignored. ^ s (x), 
s =|, | being the spin index, is the two-component electron wave function for non- interacting electrons moving to the 
right (V'ls) and to the left ("02s): 



From the Hamiltonian, we obtain the Dirac-type equation for the electron wave functions 



(2) 



( r* d 
dx 



\ 



A(x) 



and the self-consistent gap equation 



A(x) 

ih,VF~r 
dx 



Ai(x) = A(x) - A e = ^ ls (x)i,* 2 M + ^s{xWi s {x)] 



1pls(x) 



= E 



4>2s (x) 



(3) 



(4) 



where the prime on the summation symbol indicates the sum is over all occupied electronic states. 

We distinguish between the two cases with A e = and A e > 0. Let us first consider the case with A e = 0. Then 
the Hamiltonian has two degenerate ground states with A(x) = A,-(.r) = ±Ao, where Ao (> 0) is the self-consistent 
uniform Peierls gap parameter. In this paper, we focus on non-uniform solutions for the gap function. The best-known 
ones are the single soliton and the single polaron. Here we discuss the soliton and put all the algebra for the polaron 
in Appendix A. The spatial variation is set by the coherence length, a scale defined by the Fermi velocity and the 
uniform gap: 



Hvf 

~aT 



(5) 



A second dimensionless parameter k is important for polarons and is defined and discussed in Appendix A. The spatial 
variation of the order parameter for a single soliton located at x = is 



A(aO 



tanh 



The resulting (unnormalized) wave functions for unbound states with \E\ > Aq can be written as 



(G) 



ipu(x) = 4>i{x) 

1p2s(x) = (j>2{x) 



ik'.l 



, , x\ . E 

tanh — — i— ikf 

1 U Ao 



, i x \ E 

tanh — + i—. ikf 

U A 



where the quantum number k is related to the energy eigenvalue by E 2 = (hypk) 2 
explained in the next section, we introduce another equivalent set of wave functions 



A ^ 



(7) 



221. For a reason to be 



2 



-ikx 



tanh 



E 

-i—+ikf, 



fas (x) = <f> 2 (x) = e 



— ikx 



ikf, 



(8) 



which are trivially obtained from Eq. (Q) by replacing k with —k and multiplying the wave functions by i. We omitted 
the spin indices in the notations 0i, 02, 0i and 02 because the wave functions are spin- independent. There is also a 
soliton bound state located at E = and described by the wave functions 



1 



=sech 



*p2s (x) 



=sech 



(9) 
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Though there is no electronic state for < \E\ < Aq, we can still use the wave functions (0) and (||) for computing 
the subgap Green functions. 

When the extrinsic component of the gap function, A e , is non-zero, the double-degeneracy of the ground state is 
broken and there appears a higher-energy metastable state with a uniform gap. For A e > 0, the constant intrinsic 
gap Aj can be taken to be positive in the Peierls ground state and negative in the metastable state In both 

states, the uniform Peierls gap parameter is given by Ao = |A» + A e |. We caution the readers that we will use the 
same notation Ao and the definition (||), regardless of the value of A e . It turns out that the soliton solution is not 
allowed in non-degenerate cases because the equations (||) and (0) do not satisfy the gap equation (||) for non-zero 
Ao. 



III. GREEN FUNCTION AND THE DENSITY OF STATES 



In this section, we compute the retarded and advanced 2x2 matrix Green functions G + and G~ associated with 
Eq. (|J). For that purpose, we introduce two linearly- independent wave functions ip and tp that satisfy Eq. (0) in the 
interval — L < x < L and the boundary conditions that ip(—L) and ip{L) do not diverge as L becomes large, whereas 
ip(L) and ^>{—L) diverge jlB). We write the Green function in the form 



G(x,y\E) 



G n (x,y\E) G 12 (x,y\E) 
G 21 (x,y\E) G 22 (x,y\E) 



ih{x)il} 2 {y) ■ih(x)ipi(y) 
hv F {i>i$2 - i>2i>i) V ^2{x)ih{y) <tp 2 (x)i/ji(y) 
i ( ipi(x)ip 2 (y) ipiix^iiy) 



nvF^i^i - *p 2 ipi) \ 4>2(x)ip 2 (y) tp2(x)tpi(y) 



if x > y 
if x < y 



(10) 



where {^>iip 2 — ip 2 ipi) is a constant independent of x |15| , [23[|24| . We obtain G + (G ) by solving for wave functions 
with Imi? a small positive (negative) number. 

The wave functions for the single soliton case, (0i, 4> 2 ) and (0i, (f> 2 ), defined by Eqs.(7]) and (||) do not satisfy the 
necessary boundary conditions, and therefore cannot be used in the representation ([10). However, it is possible to 
construct two new solutions {ipi,ij} 2 ) and (tpi,^), which satisfy the boundary conditions, by a linear superposition 
of (01,^2) and (0i, 2 ): 



tp2 



02 



e ifeL 0i 



_ ik(x+L) 



e ifeL 2 



e- lfci 0i 

X 

tanh ( — 



e~ lkL 4> 2 



X 

tanh I — 



i— ikf 



-ik{x+L) 



X 

i aiiii ( — 



, e -ik{x+L) 



X 

tanh ( — 



e" lfcL 0i 



_ ik(x-L) 



e 02 



e ifeL 0i 



tanh 



ikL 7 

e 02 



— i— ikf 

A 



+ ie-' lk{x - L) 



tanh 



-ie 



ik(x-L) 



tanh 



i— ikf 

A 



-zk(x-L) 



tanh 



E 

i—+ik£ 



i- — h ikf 
Ao 



E 



. E 
''A~o 



ikf, 



(11) 
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First, we calculate the Green function for E > Ao. We introduce a small imaginary part E to the energy and 
replace E by E + iE. Then the quantum number k acquires an imaginary part and is replaced by k + ik such that 
EE — (hvp) 2 kk. If we choose Kvpk = — ^JE 2 — Aq < E and k have opposite signs. In order to obtain the 
retarded Green functions, we assume E > and take the L — ► oo limit first to simplify Eq. (O) using the fact that k 
is negative and e kL vanishes in this case. Finally, we take the E — > limit to obtain the simplified wave functions for 
use in calculating the retarded Green functions: 



V'i 
i>2 



a ik(x+L) 



, I X \ ■ E 

tanh — — i— ikf 

1 U Ao 



— —j P ik ( x + L ) 



4>i = ie~ lk{x - L) 



— f,—ik{x-L) 



tanh — + i—. ikf 

1 U A 

I 3l\ E 

tanh 1 — 1 — i— h ifc£ 

€/ A 

, , x\ E 

tanh — + i— h iKt 

1 U A 



(12) 



Similarly, for the advanced Green functions, we assume E < and obtain 

V-a = ie- ifc(3;+L) tanh ^ - i A + ik£ 

tanh — + i-r h 

U/ A 



V>2 
V>2 



-ik(x+L) 



ik(x-L) 



, X \ ■ E 

tanh — — i— ike 

1 U Ao 



, x \ E 

tanh — + i-T ifct 

A 



(13) 



These wave functions satisfy the required boundary conditions that tpi^(L) and ipi.2(—L) diverge as L goes to infinity, 
whereas t/'l.aC - i) and ipi^L) do not diverge. Substituting Eqs. ( |l2| ) and (|l^) into Eq. (|To|), we find 



G ± (a;,a;') 



T-2 



exp 



±iy/(E/A ) 2 -l 



Tlvf 
x{ 2 



4a/(£/A ) -1 



Ao 
Ao 
E 



1 + (7i ± V (S/A )- - 1 sgn (z - x') as 

(<T2 - 1) 



1 - t t' Q T i^J(E/A ) 2 -l \t - t' \ 

+ (t + t' Q - 2) (71 + i (t - t Q ) (7 3 



where 1 is the unit matrix and to and t' a are defined by 



I'd = tanh ( —j , t'o = tanh ^— 



(14) 



(15) 



When E < — Ao, we choose TiiTpfc = i/i? 2 — Aq |p25|] . Then E and fc have opposite signs. A calculation similar to that 
for E > Aq leads to the symmetry relationship 



Gf 1 (x,x'\E) = -Gf 2 (x,x'\-E), Gf 2 (x,x'\E)=G* 1 (x,x'\-E). 



(16) 



The subgap Green function for < \E\ < Ao will be derived in Appendix B. 

With the exact Green function in hand, it is a straightforward matter to compute the local electronic density of 
states (per spin) p s (x,E): 
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p s {x,E) = Im Tr G + (x, x\E) 

7T 



1 



1 



= \ T^hvp y/E 2 - A 2 




\E\- 



2\E\ cosh 2 (x/0 



if \E\ > A 
if < \E\ < A 



(17) 



The spatially-averaged density of states p s (E) is given by 



Ps(E) = ^ y dx p s (x, £7) 



Po(E) + Sp.(E) if|£|>A 
if < 1^1 < A 



(18) 



where the density of states in the Peierls ground state, po(E), and the correction to the density of states in the 
presence of a single soliton, Sp s (E) , are 



Po(E) 



1 



\E\ 



Sp s (E) = - 



ttTlvf y/E 2 - Ag 
1 



A 2 



2irhv F \E\ V / E 2 -A 2 L 



— 1 a uh ( — 



(19) 



In the polaron case, we derive the local electronic density of states p p (x,E) for \E\ > Ao from the Green functions 
derived in Appendix A: 



p p (x,E) 



1 



\E\ 



■xhvF yjE 2 - A 2 



X 1 



2{E 2 -E 2 ) 



1 



cosh (k(x + xq)/£) cosh (k(x — xq) /£) 



(20) 



where Eq = AoVl — k 2 and Xq is defined by Eq. (Al). The density of states when < |J5| < Ao and E ^ ±Eq is 
zero. The spatially-averaged density of states, p p (E), for \E\ > A is given by 

p p (E)= P0 (E) + 6 Pp (E), 

1 \E\ kA 2 



Spp(E) 



ttHvf y/E 2 - A 2 (E 2 - E 2 ) 



2L 



tanh 



k (L + x ) 



tanh 



k(L-x ) 



(21) 



IV. SOLITON, POLARON AND MULTIPOLARON ENERGIES 



A. Self-consistency in the uniform case 

In this subsection, we reformulate the self-consistency equation (^J) in the uniform case with a constant gap function 
in terms of the Green functions. First, we write the total energy (Et) per unit length in terms of the electronic density 
of states per spin p(E, Ao): 

where the cutoff energy E c , introduced to overcome the well-known difficulty of the continuum model that its energy 
spectrum is unbounded, is a function of the gap parameter Ao = |Aj + A e | and satisfies the requirement that the 
integrated density of states, 
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iV(£ c (A ),A ) 



dE p(E,A ), 



(23) 



is a constant independent of Ao. As mentioned earlier, in non-degenerate cases with A e > 0, the intrinsic gap Aj can 
be taken to be positive in the ground state and negative in the metastable state. Equivalently, the absolute value |Aj| 
satisfies |Aj| = A — A e in the ground state and | Z^.^ | = A + A e in the metastable state. Applying the self-consistency 
condition 



8 



dA 



E T (A ) 



21 







to Eq. ( |22| ) and using the relationships |2j: 

dN(E,A ) 



dE 



p(E, A ), 



dA 



N(E C {A ),A )=Q, 



dN< f A A ° ] = -l^[GUx,x\E)+GUxME)] 



we derive a new form of the self-consistency equation, 

| Ai| = —nfrvpA 



°° dE dN(E,A ) 



dA 



in a straightforward manner. Since in the uniform case 

G+ 2 (x,x\E)=G+ 1 (x,x\E) = 



hv F X / dEhn.[Gf 2 (x,x\E) + Gt 1 (x,x\E)] 



A 



2hv F ^/E 2 - A 2 ' 



we obtain 



|Aj| = A 



dE 



A 



AAq cosh" 



VE^AI 

'Ec 
Ao 



AA In 



E c + y/E* - A 2 
Ao 



(24) 



(25) 



(26) 



(27) 



(28) 



which reduces to 



E c = A cosh 



(29) 



in the degenerate case where Ao = |A,-|. One can easily see that, in the weak coupling limit, Eq. (£8|) is approximated 
by 



Aq = 2E C exp 



it 



Ae 

Ao 



(30) 



where the upper (lower) sign corresponds to the ground (metastable) state. The total energy per unit length, obtained 
from Eq. (22), is given by 



Ej 1 



A 2 



A 2 



cosh 1 



2L itTivf^ nhvF 
which, using the gap equation (Eq), can be rewritten as 



A 



■ sinh 



Et _ Af 
2L TrhvpX irhvp 



A 2 



AA 2 



— sinh 



2 cosh 



2|A t | 
AA 



Ec 
Ao 



(31) 



(32) 



G 



B. Soliton energy 



We compute the soliton excitation energy E 8 in the degenerate case using the density of states expression obtained 
in Sec. II] . Regardless of the number of electrons occupying the midgap state with E — 0, E s is given by 



E s = I dx A 2 [tanh 2 (z/0 - l] 

TThVpA J_ L 



4L 



r 

J An 



dE Ep s (E) 



dE Ep {E) 



(33) 



where E* is the energy cutoff in the presence of a soliton defect, while E® (= Ao cosh(l/A)) is the cutoff in the Peierls 
ground state. E s c and E® are related by the condition that the total number of electronic states in both cases has to 
be the same: 



\ + 2L 



dE p s {E) = 2L dE p (E) 



(34) 



where the number 1/2 on the left-hand side accounts for the midgap state (per spin). Substituting Eqs. (18) and ([H 
into Eq. (t}3) and using the fact that E% differs from E® by a term of order 1/L, we derive 



A A 



— tanh | — | cot. 
2L V A. 



sinh 



O 



L 2 



(35) 



Using this in Eq. 



we find, in the L — > oo limit, 



E s = — Ao cosh 

7T 









(£ 


cot -1 


sinhQ) 









(36) 



which reduces to the well-known value of the soliton excitation energy in the weak coupling limit, 2A /7r, as A goes 
to zero and increases monotonically to Ao as A grows to infinity. In Fig. 1, we plot the energy of a single soliton- 
anti-soliton pair, 2E S , together with the polaron energy in the degenerate case to be obtained in the next subsection 
(Eq. as a function of A. 

C. Polaron and multipolaron energies 

The polaron excitation energy E p is obtained in a manner similar to that in the previous subsection. Defining the 
electron occupation numbers of the subgap states at E = Eq and — E as n + and n_ (n + , n_ = 0, 1, 2), we have 



E p = —!— - / dx A? 

■nnvpA 



1 t-) 



- 1 



72-)Eo - 4L 



E l r E ° 

dE Ep p (E) - I dE Ep (E) 



An 



(37) 



where the constant Aj (= Ao — A e > 0) satisfying Eq. (|2^) is the intrinsic gap in the ground state and E% is the energy 
cutoff in the pre sen ce of a polaron defect. We emphasize that the coherence length £ appearing in the definitions of 
t+ and <_ (Eq. (A3)) is defined by Twf/Aq, not by Hvf/A-i- The equation (p3) is replaced by 



1 + 2L I C dE p p {E) = 2L [ " dE p {E). 

Ja Ja 



(38) 



Substituting Eq. ( |2l| ) into the above equation, we find 

A,; 



EP c E c , 

-r^- = — -7- tanh 
A A L VAA 



cot 



- sinh [ — r- 
k V AA 



(39) 
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The polaron energy E p to order 1 follows from this and Eq. (37) and is equal to 



4 / Aj 

E p = (n + - n_ + 2)E Q + ~A cosh — — 

7T V AA 



cot 



— sinh ( — r- 
K 1 AA 



— E tan 1 



E lAAo 



-A r (tantT 1 « - re) 



(40) 



where T = A e /AAo- The self-consistency condition is equivalent to the condition that E p has to be minimized with 
respect to n for given values of n + and Applying a weaker condition that E p takes an extremal value, that is, 
dE p /dK = to Eq. <M), we obtain 



tan 



tan# coth 



Ttantf = - (n+ - n_ + 2) , < < -. 

4 — _ 2 ' 



(41) 



where 9 is defined by k = sin 9 and i?o = Aocosf?. It turns out that, in the degenerate case with T = 0, stable 
solutions which satisfy Eq. ( pi] ) and minimize Eq. j40| ) exist only for the electron polaron state (n + = 1, n_ = 2), 
with the effective charge Q = 2 — n + — ra_ = — 1 and the spin S = 1/2, and the hole polaron state (n + — 0, n_ = 1), 
with Q = +1 and S = 1/2. In both cases, 9 = tan _1 [tanh(l/A)] and 



tanh(l/A) 



E n = 



A 



Y^tanh^l/A) + 1 ^tanh^l/A) + 1 



E p = — A cosh(l/A) cot" 

7T 



ycosh 2 (l/A) + sinh 2 (l/A) 



(42) 



which reduce to k = 1/V%, Eq — and E p — 2y/2Ao/n in the weak coupling limit. The behavior of E p as a 

function of A in the degenerate case is shown in Fig. 1. 

In the non-degenerate case with T > 0, stable polaron and multipolaron (bipolaron, tripolaron, quadripolaron) 
solutions exist for all possible combinations of n + and n_ except (n + = 0, n_ = 2). There are three bipolaron solutions 
with the effective occupation number N = n + — n_ + 2 = 2. Both the electron-electron bipolaron (n + = n_ = 2, 
Q = —2) and the hole-hole bipolaron (n + = n_ = 0, Q = +2) have spin zero, while the electron-hole bipolaron 
(n + = ri- = 1, Q = 0), which can also be called an exciton, exists in either spin singlet (S = 0) or triplet (S = 1) 
forms. More details about the interpretation of multipolaron solutions can be found in Refs. [0J1. Explicit values 
of k, En and E v for given values of A and T can be obtained only numerically. As a specific example of applying 
Eqs. (|40|) and (|4J, we consider the case of c«s-polyacetylene. There is a considerable uncertainty in the precise values 
of A and T (= A e /AAo). In order to make a rough estimate of the excitation energy, we use a set of the approximate 
values close to those assumed in plL which are A = 0.4, A e — 0.06 eV, Ao = 1 eV, and therefore, T = 0.15. The 
results are listed in Table I. 



V. CONCLUSION 



In the present paper, we have reformulated the theory of solitons, polarons and multipolarons in both degenerate 
and non-degenerate conducting polymers. Especially, we have developed a simple method for calculating the Green 
function and the density of states in the presence of a soliton or polaron defect and computed the soliton, polaron and 
multipolaron excitation energies and the self-consistent gap functions exactly for an arbitrary value of the electron- 
phonon coupling constant (See Eqs. ( ^6| , ^p[^ ) ) . Our method can be generalized in a straightforward way to more 
complicated situations. In a previous work, Kim and Wilkins have devised an efficient numerical method for calculating 
the exact disorder-averaged Green function in the presence of a soliton or polaron defect and disorder fl5|| . It appears 
that, by combining this method with the results of the present work, one can study the disordered soliton or polaron 
problem in an exact and self-consistent manner. Work in this direction will be presented in a separate publication. 
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APPENDIX A: POLARON GREEN FUNCTION 



The spatial variation of the order parameter for a single polaron located at x — is given by 

A(s) 
A 



xo 1 . 
-r = i— In 

£ 4k 



I - ; . s tanh 

1 + K 



k(x + x ) 


— tanh 


/c(x — x ) 


} 


[ £ 




[ t . 





1 - K 



< K < 1. 



(Al) 



The dimcnsionlcss parameter k appearing in this equation determines the shape of the polaron, which changes from 
an extremely shallow well centered at x = for < n <C 1 to a well-separated pair of soliton and anti-soliton located 
at xq and — xo, respectively, as k — > 1. This parameter cannot take an arbitrary value, since the self-consistent gap 
equation is not satisfied for every value of n. When A e = and in the weak electron-phonon coupling limit (that is, 
the A — > limit), it turns out that k has the unique value of 1/V2, as will be shown in Sec. IV C . 

In case of the single polaron configuration, the exact unnormalized wave functions for unbound states with |J5| > Ao 
can be written as [E2| 



E 
Ao 
E 

Aq 



- 1 + 7(1 + t)i+-<y(l- t)t_ 
+ 1-7(1 -i)t+ + <5(1 + *)t_ 



(A2) 



where t+, 7 and S are defined by 



tanh 
1- 



k(x + Xq) 



ifikpk 
E + A 



tanh 



k(x — Xq) 



ihkpk 
E + A 



(A3) 



Similarly to the soliton case, we introduce an equivalent set of wave functions 



ipu(x) = 4>i(x) 



^2s{x) = 4>2{x) 



-ikx 



-ikx 



E 

k£ - -r- + 1 - 5(1 + i)t+ + 7(1 - i)t. 

Aq 



E 



1 + 6(1 - i)t + - 7(1 + i)t. 



(A4) 



which are obtained from Eq. (A2) by replacing k with —k and multiplying the wave functions by —1. In addition, 
there are two subgap polaron bound states located at E = ±Eq = ±Aovl — k 2 , the wave functio ns f or which we do 
not write down here explicitly. As in the soliton case, we can also use the wave functions ( |A2|) and (A4) for computing 
the subgap Green functions when < \E\ < A and E ^ ±E . 

We superpose two independent solutions in the polaron case, (A2) and (A4), to make two new wave functions 
(ipi,ip2) and (?/>i, ^2) satisfying the necessary boundary conditions. When E > Aq, following the same procedure as 
in Sec. Ill, we obtain the simplified wave functions 



>h 

V'l 

V'2 



_ e ik(x+L) 
_ ^ik(x+L) 
— g-ifc(>-L) 



E 



E 



1 + 7(1 + »)*+- 5(1 -i)t- 



k£ - -r- + 1 - 7(1 - i)t+ + 5(1 + i)t. 
Ao 

k£- -^- + 1-6(1 + i)t++ 7(1 -i)t. 

k£+ — -1 + 5(1- i)t+ - 7(1 + i)t- 
^o 



(A5) 



with hvpk = —^E 2 — Aq |g5[ for computing the retarded Green functions and 



9 



fa = e -iH*+L) 

fa = e - lk(x+L) 
fa = e M*-L) 

fa = e ^- L ) 



k£-^- + l-5{\ + i)t+ + 7 (1 - i)t- 

k£ + — - 1 + 7 (1 + i)t + - 5(1 - i)t_ 
£;£ - X" + 1 - 7(1 - »)*+ + 6 ( 1 + *)*- 



(A6) 



for the advanced Green functions. The Green functions obtained by substituting Eqs. ( A5) and (A6) into Eq. (|l 
are 



G ± (a;,a;')= 



exp 



±iyJ(E/A f-l \x-x'\/£ 



flVl- 



aJ(e/^Y-\ 



E 



< 2 ( — 1 + g x ± V (^/A ) 4 - 1 sgn (s - a;') ct 3 



* 2 Ag 
E*-E* 



1 - t_tl T *V (-E/Ao)' - 1 |*_ - (1 - <7 2 



E 
A^ 

1 - M'+ T i\[{E/&o? - 1 |*+ - J (1 + <r a ) 
(l - (£/A ) 2 ) (t_ + t'_ - i+ - t' + ) /k + 2 - - t+t'_ 
T i\f(E/A f -l(t+t'_ - t-t' + + (t- - t'_ +t+- t' + ) /K^sgn(a; - x') 

Ti\l{E/Aof "if 2 - *-*+ - M- + (*- + *- -*+ - *+) /«^sgn(a; - x') 
+ (l - (^/A ) 2 ) (t- - f'_ + i+ - 4) l K + M- - *-*+ 



(A7) 



The Green functions for E < —Aq are obtained using the symmetry relationship Eq. (jig), which is valid in the polaron 
case, too. 

When < \E\ < A and E ^ ±E , we obtain 



fa = e q{x+L) 
fa = e q{ - x+V > 
fa = e - q{x - L) 
fa = e~ q{x - L) 



E 

+ T- - 1 + 7(1 + i)t+ - 5(1 - 

- -r- + 1 - 7(1 - + <5(1 + »)*- 

- -r- + 1 - 5(1 + + 7(1 - *)t- 

+ -r- - 1 + 5(1 - l)t+ - 7 (1 + i)t- 
^0 



(A8) 



where q = ik = y 1 — (E/Aq) /£ (See Appendix B) and 7 and 5 are given by Eq. JA3Q with zfc replaced by g. 
Substituting Eq. @) into Eq. ©, we obtain 



exp 



G ± (x.x'y 



■y/l-{E/A f\x-x'\/t 



flVj- 



Adl-iE/AoY 
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x 1 2 ^1 + a x + 2 ^/l-(£/A ) 2 sgn {x - x') a^j 



* 2 Ag 
Et-E* 



E 



1 - t-i 7 . + V 1 - (E/A y\t- - t'_\/K (1 - CT2 ) 



+ l - £+t^ + V 1 - ( E / A o) \t+ - t'+\/K (i + as) 



(l - (£/A ) 2 ) (*_ + *L - f+ - t+) /« + 2 - t_t' + - t+t'_ 

1 - {E/A a ) 2 U+t'_ - t_t' + + (t_ - + 1+ - sgn (a; - x') 

yi - (E/A a ) 2 (^2 - - t + t'_ + (i_ + i'_ -t+- t' + ) /«^sgn(a; - x) 

+ (l - (B/A ) 2 ) (t- - +t+- t' + ) /k + t + t'_ - t_t' + 
where t' + and t'_ are denned similarly to t+ and t_ except that x is replaced by x'. 



(A9) 



APPENDIX B: SUBGAP GREEN FUNCTION IN THE SOLITON CASE 



Here we calculate the Green function for < \E\ < Ao. In this case, we have kt; = Hy/l — (Ej Ao) 2 = ±ig£. We 
choose k = —iq and take the L — > oo limit to simplify Eq. (11) using the fact that g is positive and e~ qL vanishes: 



tanl>(!)-iiL-tf 



V^ 2 = -ie 9(x+i) 
V3 2 = e -9( a; - L ) 



tanh 
tanh 



tanh | - J + i-^- + r/i 









(!) 


A 




E 




(!) 





The Green functions are obtained in a straightforward manner by substituting Eq. (Bl) directly into Eq. 
find that the retarded and advanced Green functions are the same in the present case and given by 



G ± (x,x')= 



1 



exp 



-yJl-(E/A ) 2 \x-x'\/Z 



xl2 



hvF 

' E 

Ao 

Ao 

E 



AJl - (E/A ) 



1 + a x + id 1 - (E/A ) sgn (x - x') a 3 



l-t ^ + Jl-(E/A y\t -t' 



(<T2 - 1) 



+ (t + t'o - 2) a x + i (t - t' ) a 3 



(Bl) 
. We 



(B2) 
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FIG. 1. The energy of a soliton-anti-soliton pair, 2E S and the polaron energy in the degenerate case, E v vs the electron-phonon 
coupling constant, A. 



TABLE I. Polaron and multipolaron states in cis-polyacetylene when A = 0.4 and F = 0.15. N = n+ — n- + 2 is the effective 
occupation number, Q = 2 — n+ — n_ is the effective charge, 6 is the angle (in radians) defined by Eq. (^) and E p is the 
excitation energy. 



n+ 


n_ 


N 


Q 


e 


E P /Ao 


Interpretation 


1 


2 


1 


-l 


0.66 


0.93 


electron (e) polaron 





1 


1 


+1 


0.66 


0.93 


hole (h) polaron 


2 


2 


2 


-2 


1.19 


1.51 


e-e bipolaron 


1 


1 


2 





1.19 


1.51 


e-h bipolaron (or exciton) 








2 


+2 


1.19 


1.51 


h-h bipolaron 


2 


1 


3 


-1 


1.41 


1.75 


e-e-h tripolaron 


1 





3 


+1 


1.41 


1.75 


e-h-h tripolaron 


2 





4 





1.48 


1.87 


e-e-h-h quadripolaron 
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